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Abstract 

We provide evidence that the small-scale statistics of the accelera- 
tion of a test particle in high-Reynolds number Lagrangian turbulence 
is correctly described by Tsallis statistics with entropic index q = |. 
We present theoretical arguments why Tsallis statistics can naturally 
arise in Lagrangian turbulence and why at the smallest scales q = | 
is relevant. A generalized Heisenberg-Yaglom formula is derived from 
the nonextensive model. 
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Recently, methods from nonextensive statistical mechanics [[I], |3| have 
been successfully applied to fully developed turbulent flows [|J |], |, |7j. 
As a driven noneqilibrium system a turbulent flow cannot extremize the 
Boltzmann-Gibbs entropy — it is obvious that ordinary statistical mechanics 
fails to provide a correct description of turbulence. But there is experimental 
and theoretical evidence that the statistics of velocity differences is well de- 
scribed if one assumes that the flow extremizes the 'second best' information 
measures available. These are the Tsallis entropies [p]] defined by 
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q 7^ 1 is the entropic index and Pi are the probabilities associated with the 
microstates of the physical system. The Tsallis entropies, closely related 
to (but different from) the Renyi information measures ||, |9|, are convex 



and distinguished by generalized Khinchin axioms [10]. For q — > 1 they 



reduce to the ordinary Boltzmann-Gibbs entropy. Their importance has been 



demonstrated in numerous recent papers (see |TT| for a detailed listing). 

Whereas previous papers on turbulence and Tsallis statistics mainly dealt 
with the inertial range j|, §|, in this Letter we will concentrate on the small- 
scale characteristics of fully developed turbulent flows. Here one is still far 
away from a complete theory, though many empirical facts of the small-scale 



statistics are well known and have been experimentally verified (see e.g. |T2 
for a review). Most turbulence measurements have been conducted in the 
Eulerian frame, i.e., by measuring the spatial fluctuations of the velocities 
using the Taylor hypothesis. Recently, experimental progress has been made 
in investigating the Lagrangian properties of fully developed turbulence, by 
tracking test particles that are advected by the flow [13|, [14]]. Examples 
of measured histograms of the acceleration of a test particle as measured by 
Bodenschatz et al. are shown in Fig. 1. The distribution has been rescaled to 
variance 1. The acceleration has been extracted by parabolic fits over 0.75t v , 
where = {v/e)^ is the Kolmogorov time, v is the kinematic viscosity and 
e is the average energy dissipation rate. The figure also shows as a solid line 
the normalized probability density function 

p(x) = -(1 + 2x 2 + x 4 )- 1 = -(1 + x 2 )- 2 , (2) 
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which has variance 1. Apparently, the experimental data and the above 
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distribution function agree quite well. We will now present theoretical argu- 
ments why this type of distribution is relevant. 

The acceleration a of the Lagrangian test particle is a strongly fluctuating 
random variable. It can be regarded as a velocity difference on a small time 
scale r, i.e. a « (v(t) — v(t + r))/r := u/t, where r := kt v is of the order 
of the Kolmogorov time. It is not clear whether the limit r — > exists 
in a mathematically rigorous way. For example, already for an ordinary 
Brownian particle described by the Ornstein Uhlenbeck process the 
velocity exists but the acceleration is singular. Hence it seems to make sense 
to statistically describe the movement of the Lagrangian test particle using 
a small effective finite time scale r. One can now develop the formalism of 
nonextensive statistical mechanics for the fluctuating temporal changes in 
velocity. Following the ideas of @] we define formal energy levels E(u) by 
the kinetic energy of the velocity differences 

E{u) : = \a 2 r 2 « l(v(t) - v(t + r)f = \u\ (3) 

Moreover, a formal temperature Pq 1 is introduced as 

p- 1 := 6T = eKT v = Ke l,2 u 1 ' 2 (4) 

The multiplication with a time scale r is necessary for dimensionality reasons, 
since e has dimension length 2 /time 3 . Extremizing the Tsallis entropies one 
obtains the following generalized version of a canonical distribution 0] 

K«) = ^ (1 + ^(9- " 9 - (5) 

For q — > 1 the above probability density reduces to the ordinary Boltzmann 
factor p(u) ~ e -2^ u . The value of the entropic index q > 1 depends on 



the Reynolds number and the scale (see |19| for precision measurements in 



Eulerian turbulence). The normalization constant Z q is given by 

For q — I and if rescaled to variance 1, the distribution (|5|) is identical with 
the distribution (Q), which apparently is in very good agreement with the 
experimental data. 
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All moments of the generalized canonical distribution (|5]) can be evaluated 
analytically. In particular, one obtains for the second moment 
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<«>> = y__ .■*.)*.= m 

This yields for the second moment of the acceleration a = u/r 

(a 2 ) = n \ 2 = e^u- 1 / 2 -^— (8) 
X 1 (3 T 2 5-3q K5-3q K J 

Thus we obtain from the nonextensive model the Heisenberg-Yaglom relation 
(a 2 ) = aoe 3 / 2 ^ 1 / 2 , identifying the constant ao with 
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Gaussian statistics (q = 1) and k « 1 would imply a m 1. On the other hand, 
the true turbulent small scale statistics is much better described by q = 3/2, 
which yields ao ~ 4, in agreement with direct numerical simulations for large 
Reynolds numbers |Tj|. The precise numerical value of ao also depends on 
the ratio k = t/t v which enters into the formal thermodynamic description 
via eq. (f|). Bodenschatz et al. measure distributions with q w 1.5 and 



ao ~ 5.3 for large R\, which corresponds to k ~ 0.75. Generally, a measured 
Reynolds number dependence of a (as presented in |H|) can be translated 



into a measurement of the entropic index q in Lagrangian turbulence, by 
solving eq. @ for q. 

Let us now argue on theoretical grounds a) why Tsallis statistics can 
naturally arise in Lagrangian turbulent flows and b) why the entropic index 
is q — | at the smallest scales. 

Generally, Tsallis statistics with q > 1 can arise from ordinary statistical 
mechanics (with ordinary Boltzmann factors e~^ E ^) if one assumes that the 
formal temperature (3~ l is locally fluctuating (see [17] for similar ideas). In 
our application to a turbulent flow, is identified with the product er of 
local energy dissipation and the typical time scale r during which energy is 
transferred. Both quantities can fluctuate. is a formal variance param- 
eter describing the fluctuating environment of the Lagrangian test particle, 
measured relative to the movement of the particle. Using the integral repre- 
sention of the gamma function 



-H z - l dt (10) 
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and substituting 



(12) 



q-l 



one may write 
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1 + (q - l)(3 Q E{u))-^ = / e-WfWdp (13) 
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with 



/( ^F^I(^P A M-(^kl (14) 

being the x 2 distribution. The physical interpretation of eq. fll3|) is that 
due to fluctuations of (3 with probability density f{(3) the Boltzmann factor 
e -pE{u) Q £ or dinary statistical mechanics has to be replaced by the gener- 
alized Boltzmann factor (1 + (q — 1)[3qE{u)Y~ of nonextensive statistical 
mechanics. The Tsallis distribution with fixed variance parameter /3q effec- 
tively arises by integrating over all possible fluctuating variance parameters 
(3. This illustrates why the nonextensive formalism can be relevant to non- 
equilibrium systems (formally described by a fluctuating 0) if there is a 
quasi- stationary state in probability space. 

The x 2 distribution is well known to occur in many very common cir- 
cumstances (see e.g. ||18||). For example, if one has n Gaussian random vari- 
ables Xi then the sum of the squares of their deviations from their mean 
X 2 ■= i Er=i( X * - X? is X 2 distributed with ^ = n - 1. Or, if T = 
obeys a linear Langevin equation with a constant source term and a damping 
that fluctuates stochastically, one also obtains a x 2 distribution for (3 Jl7| . 

The average of f3 is given by (3 = J" °° (3f(f3)d{3 = (3$ and the variance by 

2 

(3 2 — (3 — (q — 1)/3q. From this one obtains a physical interpretation of the 
entropic index q in terms of the variance of (3, namely 

q=L (15) 

P 
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If there are no fluctuations of (3 at all, as in ordinary statistical mechanics, 
eq. fljjD just reduces to q = 1, as expected. 

In a turbulent flow, the variance parameter (3 surrounding the Lagrangian 
test particle fluctuates, and hence Tsallis statistics can naturally arise in this 
context. Since the fluctuations of j3 become smaller if the volume r 3 over 
which the energy dissipation is averaged increases, q must be a montonously 
decreasing function of the scale r. In fact, at largest scales r one observes 
approximately Gaussian behaviour (g« 1), in the inertial range q = 1.1. ...1.2 
gives good fits of the experimental data of Eulerian turbulence @, [19]], and 



Fig. 1 indicates that q ~ 1.5 at the smallest time scales of Lagrangian tur- 
bulence. 

Let us now provide a theoretical argument why q — | at the smallest 
scales. The observation is that for large \u\ the Tsallis distributions (|5|) (also 
called student or t- distributions in the statistics textbooks ) decay as 
where w = -^j. Hence only moments (|u| m ) with m < w — 1 exist. If q — | 
at the smallest scale, this means w = 4 and hence the third moment would 
just cease to exist at the smallest scale. If q is precisely | the third moment 
is logarithmically divergent, if q = | — + it just exists. Since generally the 
third moment is the most important moment in turbulence, related to average 
energy dissipation, the existence of this moment is necessary for turbulence 
to make sense. Since q is monotonically decreasing with scale r, one ends up 
with the largest allowed value of q at the smallest possible scale. This is just 
q= |-0+. 

There is further experimental evidence for the above conjecture on the 
small scale statistics. In [l9j systematic measurements of the exponent w(r) 
were performed for Eulerian turbulence (for a turbulent Taylor Couette flow). 
The measurements were performed at distances r much larger than the Kol- 
mogorov scale rj. Over a large range of scales r the measured exponents w(r) 
were very well fitted by a power law of the form 



w (r) = 4 I -J (16) 

with 5 ~ 0.3. Extrapolating this down to the Kolmogorov scale rj one obtains 
w = 4 at r = rj, which again supports our hypothesis. Although for very 
small r deviations of this power law (pointing towards smaller values of q) 
were observed in |L9], these deviations can be explained by the disturbing 
effects of noise, which naturally shifts the entropic index to lower q values, 
since Gaussian white noise implies q = 1. 
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Our conjecture on the small scale statistics is also consistent with a large 
amount of other experimental data. In |12 the Reynolds number dependence 
of the third and fourth moment of u at the smallest possible scales was 
analysed. By averaging the data of many experiments it emerged that the 

1 /3 

4th moment increases roughly like ~ i? A , whereas the 3rd moment stays 
almost constant or increases much less rapidly with R\ than the 4th moment. 
This means that at the smallest scale the 4th moment is expected to diverge 
for R\ — > oo and the 3rd moment may either just exist or may weakly diverge. 
All these experimentally observed features are correctly reproduced by the 
Tsallis distribution (|5|) with <? = §. 
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Figure captions 

Fig. 1 Experimentally measured probability density of the acceleration of a 
test particle in Lagrangian turbulence for Reynolds number R\ = 200, 690, 970, 
respectively, and comparison with the distribution 
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